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ABSTRACT 


The regulation of glucose concentration in the body has triggered 
more often the stimulant approach. The hourly operation of insulin- 
glucose-dependent bodily functions has left many type-I diabetic 
patients with a lot of pain. Doctors have proposed various measures 
such as oral ingestion of glucose to trigger insulin production to 
increase glucose uptake in the body. Such measures have not been 
mathematically modeled. In the current, a mathematical model 
formulating the concentration of glucose and insulin where extra 
glucose intake leads to extra stimulation of insulin production has 
been presented. The mathematical rigor such as the well-posedness of 
the system indicated that the proposed model is applicable in real-life 
applications. Numerical simulation indicates that both glucose and 
insulin concentration continuously increases, a sign of a need to find 
an optimal control strategy in order to manage type-I diabetes. A 
sensitivity analysis showed a constant rate of insulin-independent 
glucose disappearance is the most sensitive in the model. Optimal 
control investigation indicated that increasing the constant rate of 
insulin-independent glucose disappearance is effective after 15th 
hour. Future studies need to consider including other parameters, 
such as glucose accelerated loss due to exercise in the system. 


KEYWORDS: Oral 75g glucose intake, Glucose concentration, 
Insulin Concentration, Insulin Stimulant 


1. INTRODUCTION 

Regular human insulin production in the body begins 
to function after 30-60 min, and the peak 
concentration is witnessed after 2-4 hours [1]. The 
adverse effects are felt between 5 to 22 hours. 
However, there is no clarity on whether extra glucose 1 
usually taken by the patients effectively manages 
these effects. 


Usually, 75g of glucose is taken orally to stimulate ) 
insulin production to regulate glucose concentration 
[2]. Absurdity often occurs when the patients take 


Contribution 

The proposed study has two major contributions 
The glucose intake to stimulate insulin production 
is evaluated to ascertain its effects on the overall 
glucose concentration in the body. 


How to cite this paper: Irene Melly 
Chebet | Rotich Titus | Bitok Jacob "A 
Model of Insulin-Glucose Concentration 
for Oral Ingestion of 75g of Glucose for 


Type-I Diabetes Management" 
Published in 
International 


Journal of Trend in 
Scientific Research 
and Development 
(ijtsrd), ISSN: 
2456-6470, 

Volume-6 | Issue-7, 
December 2022, pp.363-369, URL: 
www.ijtsrd.com/papers/ijtsrd52338.pdf 


Copyright © 2022 by author (s) and 
International Journal of Trend in 
Scientific Research and Development 


Journal. This is an 

Open Access article 
distributed under the cas 
terms of the Creative Commons 


Attribution License (CC BY 4.0) 


(http://creativecommons.org/licenses/by/4.0) 


the research area [3, 4], leaving a hole in the 
mathematical perspective. 


Establishing an optimal control of the system to 
regulate the glucose concentration in the body in 
the presence of extra glucose and insulin 


extra glucose but still experience the side effect of 
excess glucose in the body, which often leads to type- 
I diabetes. 


Mathematically, the oral ingestion of glucose to 
stimulate insulin production for glucose uptake by the 
body has not been studied. The biological 
interpretation of these phenomena has not been 
widely received with open _ understanding. 
Biologically, inter-organ dependence has dominated 


production. 


The rest of the paper is organized as follows: Section 
2 presents the model formulated and governing 
equations. Section 3 the wellposedness of the system, 
investigating the system application in real life. 
Section 4 presents the equilibrium of the system, such 
as disease-free equilibrium and endemic equilibrium, 
and basic reproductive number. Section 5 presents the 
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estimation of parameter values, sensitivity analysis, 
and numerical simulation of glucose concentration 
and insulin in the body. Section 6 presents a 
formulation of the optimal control of the system. 
Section 7 presents the conclusion of the study. 


2. Proposed Study 

2.1. The model formulation 

The model presented in Figure 1 shows the 
interaction between glucose & and insulin J 
concentration in the body. Glucose is increased in the 
body at a constant rate Iz due to body metabolism 
such as food conversion. It is also increased in the 
body at a constant rate “PW; due to 75g of glucose 
intake. It is depleted at a constant rate #2. Insulin 
concentration in the body is increased due to oral 
ingestion of glucose of 75g at rates #5. The insulin is 
decrease at rate and #2 due to a constant rate insulin 
degradation and a constant rate insulin-dependent 
glucose disappearance 13. 


IT. 
32 GD 
| ce 
Pe 72 + BS 
Pa z P2 - 
¢ FT 2 


Figure 1: A model formulating the interaction 
between glucose and insulin. 


2.2. The governing equations 


ac 

a My + Wz — (og +92IG (1) 
dl 

a 2 — Pal — 71246 (2) 


with @ > 0,2 > 0. 


3. Well-posedness 

3.1. Existence of the solution 

We show that the model (1) - (2) solution is unique, 
non-negative, and bounded, making it applicable in 
real life. Let X(t) = CG(e),J(¢)) and f:X¥ =G' such 
that f = (f,, fs, where 

A(X) = Mz + Pz — (92 +21) (3) 

f,(X) = Bz — Pal — NIG (4) (5) 

So that equations (1) - (2) can be written as 


X= F(X(t)):X(0) = (Go,lo) (6) 


Theorem 3.1 Suppose f(X(t)) is given by (6), with its 
initial condition X(Q) = G +I > 0, is non-negative. 
System (1)-(2) solution is unique, non-negative, and 
bounded. 


Proof. if f; are continuous functions and oh. 
a 


1 =i,7 =2 exists then f(X(t)) is locally Lipschitz 
continuous. (0) = G+ I > 0. Therefore, at least one 


compartment is non-empty. This means that a unique 
solution X(t) of the system exists, defined in some 


time interval containing t = 0 [5]. If we let tg be the 
smallest ¢ such that G(tg) = 0 or I{tg) = 0. We use 
continuity of G(t), and I(t), 3 t* > tg such that if 
G(t;)=0, then from (1) to show _ that 
<= =I, + ¥; > 0, Wt € [tp,t*]. This means that G is 
an increasing function on interval [t9,¢*], which in 
turn means that G(f) = OWt € [tp,t*]. We use the 
same approach to show that, if J{t).) = 0, then from 
(1) we get that == @ 20 — I(t) > OVE [tp,t*] 
Thus, the solution to the system is non-negative. We 
use the dissipativity condition of theorem 2.3.6 in [6] 


to show that the solution of the system (1)-(2) exist 
globally as follows; 


f(x). X= (Cf fo) 
= Gf, +If2 (7) 
= GU, + ¥2-—G(a2 + Iy2))+ 
I(@ — Mp2 + Gyz)) 
Therefore, a unique solution ¥(#)4'vt = 0. Hence the 
solution is bounded. 
4. Stability Analysis 


4.1. Endemic Equilibrium 
The equilibrium state is defined as when there is no 


change in the system, = 0. Endemic equilibrium is 
ad ‘ . 
when = 0, but no state variable is zero, that is, 


I! =Qand G #0. We derive the endemic equilibrium 
as follows; 


We use (1) to get; 
ed 


Tits 
ct = —. 8) 
eth ( 
We also use (2) to get; 
yt, 
gt = Teer (9) 
Wat 
We use (8) and (9) to get. 
a Tens #anst¥stateces-H | 
>= (10) 
272A: 
K, = (@3n}, — 221 ,n3 — 26 ,'¥2n} + 262n2 020, + Mins + 
211, ¥n2 + 211.9, p,0, + ¥2ni + 2¥.N,p20, + piaz) 
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4.2. Disease Free Equilibrium 
The disease-free equilibrium, “ = 0 but Gy = O and 


Ij = 0.We follow algebraic oe to show 


that Gj = es IG = 2. Hint we assume the 
og tnet: 


4.3. The Basic Reproduction Number 
We follow the works of [5, 6] to compute Rg as 


follows: 


Fa (ria) Y= Gu) 


We compute the jacobian matrix to attain 


Nol ), ie 0 ) 
F= ( ¥ = i, 12 
H2f mz2G/" 0 pz oe 
We attain F at the DFE as; 
1, Y 
N2P2 Ne Sette 
Fore = *naFa | (13) 
DFE ~ Nigt+#y ) f 
i2P2 We or 


(14) 


We use (13) and (14) to compute FV~+ as follows; 
Faiz al Tats) 


(og+# ons) 
data tala?) (15) 


Gq = Paleat ana) 
We compute the eigenvalue to get Rg as; 


' Paths, (e2 . dgt Pols M16) 


2 is the constant rate of decrease of glucose 
concentration due to insulin-dependent constant rate, 
#2 is the constant rate of insulin degradation; thus, 
Pz 1s the rate of degradation of insulin- and glucose 
in the body. Iz is the glucose intake in the body, and 
@ is the constant glucose disappearance. ’; is the 
oral ingestion of glucose, and 3 is insulin increased 
due to oral ingestion of 75g of glucose. Therefore, Fg 
is the sum of the constant rate of glucose and insulin 
reduction in the body. 

4.4. Stability of DFE 

4.5. Local stability analysis of the DFE 


Theorem 4.1 Suppose Xj = 9 %2)i is a DFE 


for the system (1) and (2), then Xg is locally 
asymptotically stable if Ry «. 1and unstable if Rg > 1 
where Rg is defined by (16). 


Proof. See proof of Theorem 4.1 in [7]. 


4.6. Global stability analysis of the DFE 
Theorem 4.2 Suppose Ry = 1, the DFE of system (1) 


and(2) is globally asymptotically stable in 
={¥ =(6, €R2:G+12Cg andifRy> 1, the 
DFE is unstable. 


Proof. We use a matrix-theoretic method as suggested 
by [7]. We assume x =(G,/)? and y = (G,f). We 
consider F, !’, F and ¥ as defined in Section ??; then 


dx 
dt =F (x, ¥) —V(x,¥), 


and dynamics of infected compartments is given by 
ax 
= F-Vx — Fey), 


thus, we can obtain f{x,¥)} as follows 
f(my) = (F —V)x— F(x y) + V¥(ay¥) 


= (Ss (02 tna) Gin Gaz + Gn — Ging = 
Toa + Pn — Gin Te— 1(p,— Gna)— Gin} 
while 
crue a2( Mato) 
-ip_|{ % slot Fo7g) | 
Vit = ony Bolt Fy) 
Pz Palogteong) 
We see that F200, votsa fixyv)=0, 


r(x (Ze T,+ “=, @,)) =0 in ©. It is also clear that 


VF is pa We can therefore construct a 
Lyapunov function based on Theorem 2.1 and 2.2 as 
stated by [8]. We denote the left eigenvectors of V-+F 
corresponding to the eigenvalue of Ry by {1,172}. 


Then 
(vy, v4 )V-1F = Ro(v,, v2) 


or 
(vr, r vv “F= = 


(®22(2+2 i D eee eas 


and 


= eT o2[Tig sal 
Rofvy,v2) == (py, + ==) Wry va) (19) 


We use equations (18) and (19) to derive a possible 
solution of t 4,02. To allow work-ability, let 
1% =v = 1. Thus, by Theorem 2.1 of [8], we have 
@ = v,V¥~1txas the Lyapunov function of the model is 
given by; 
Q=v,V-1x 
ae oe 

t Pz 
We differentiate @ at DFE to get 
Qlore' = (Ro— 1) vx —v,V~1 f(x,y) 
#0 


@ IJTSRD | Unique Paper ID —- NTSRD52338 | Volume -6 | Issue—7 | November-December 2022 


Page 365 


International Journal of Trend in Scientific Research and Development @ www.ijtsrd.com eISSN: 2456-6470 


+ T ! 
This implies that x = 0 and f(x,y) =0 or y = (6,17 = (—— 2) . Therefore, (= ©.) is the only 


tists" tees" 
invariant set in 2 where Q' = 0. Thus, by LaSalle’s invariance principle [9], the DFE, 2 . ,) is globally 
‘ ta 


asymptotically stable in and Ry > 1. 


5. Numerical SImulation 
5.1. Parameter Estimation 
Table | gives a description and approximate values of the parameters. 


Parameter it I th A | Source | 
Insulin increased into the body due 
D. : 
: 75g glucose intake an ae) [10] 
Glucose increase into the body due 
yy - 
75g glucose intake ee [10] 
Ng Glucose infused into the body 1.56 1.4-1.7 [11] 
2 Constant rate insulin-dependent 3810-8 28x10-8—5.x10-8| [11] 
glucose disappearance 
P2 Constant rate insulin degradation 0.0437 0.01-0.1 [11] 
a Constant rate insulin-independent 0.0226 0.01-0.09 [11] 
glucose disappearance 


5.2. Sensitivity analysis of the parameters 

We perform sensitivity analysis to determine the most influential parameters on Fg given by (16). This will 
inform the decision on the parameter that influences the concentration of & and J. Sensitivity analysis is useful in 
showing the closeness and relatedness of the parameters. A smaller change in parameter may result in a bigger 
change in Hg. Figure 2 summarizes sensitivity analysis for the model parameters to F5. 


Sensitivity Analysis 


T T T ———————— ee ee T T 


45 4 35 3 25 2 15 4-1 0.5 ) 0.5 
«101 
Figure 2: Sensitivity analysis showing 7 is the most sensitive. Thus, the constant rate of insulin- 
independent glucose disappearance is the most influential in the £5. 


5.3. Numerical Simulations of the model system 

The numerical simulation shows the curves of the concentration of glucose and insulin. Figure 5.3 shows the 
glucose and insulin concentration curves. The curves indicate that insulin concentration increases faster than 
glucose concentration. 
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Figure 3: A plot of glucose and insulin 
concentration for a span of 48 hours. 


Figure 3a shows a continuous glucose concentration 
growth throughout the 48 hours understudy. The 
observation implies that insulin production does not 
control glucose concentration since insulin 
concentration in Figure 3b concentration increases 
until 24 hours, then the curve begins to reduce. While 
the curve for insulin concentration begins to reduce, 
the curve for glucose concentration continues to 
increase. This clearly shows a need for an optimal 
control strategy to manage glucose concentration 
levels. 


6. Optimal Control 

A sensitivity analysis of the system presented in 
Figure 2 indicates that a constant rate of insulin- 
independent glucose disappearance @2 is the most 
sensitive in controlling glucose-insulin concentration. 
Therefore, we formulate optimal control that can be 
attained by increasing or reducing the #2. However, in 
this case, the objective function is to increase #2 to 
reduce or flatten the glucose concentration at some 
point. 


If we let & be the constant rate of increase of @, then 
a set of admissible @; control is 


{@2(¢): [0, 7] = [0,1]}. (20) 


6.1. The Extended Mathematical Control 

To effectively reduce the glucose concentration, 
fz + @2 increases the overall constant rate of insulin- 
independent glucose disappearance, thus the extended 
model can be written as follows: 


= = Tl, +P2— (a2 + G2) +216 (21) 
al 
a 7 Pe Pal —12fG 


with initial G{t) > O,/(¢) > 0. A successful control 


minimizes glucose concentration. Thus, the control 
u = (@;) is considered optimal if it minimizes the 
objective function defined as 


r 
I =| (Ayg[Mz — Fz — ((02 + 0g) + 21)G] + Ag2@,)at 
? (22 


Where Ay2,A22 are the unit cost of increasing the 
insulin in the body per unit of time. 


6.2. Minimization of the objective function 

We employ Pontryagin’s maximum principle and 
optimality conditions to obtain the solution to the 
optimal control problem (23) subject to (24). 


oct 
ne (23) 


Jy (Ar [a +B 2 — (oz + 82) + 2G] + ApOa)dt 


Subject to 
= = [fz + Fy - ((u2 + i) + m21)G G(0),1(0) > 0 (24) 
_ a ®; — pol —nalG LN = 8 


A numeral graphical solution to the problem is 
presented in Figure 4. 


St) is without optimal control, G,(t) is with optimal control 
18 5 : : 1 : x 7 


16} [os-- Ga) - | 


Concentration 
% 


Time in Hours 
[a] Glucose Concentration 
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I(t) is without optimal control, I(t) is with optimal control 
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[b] Insulin Concentration 
Figure 4: Glucose and insulin concentration 
before G(t), 2(f) and Gg(t),/(f) after optimal 


control. 


Figure 4 shows the effectiveness of the control 
parameter used. A comparison of Figure 6.2a 
indicates that glucose concentration increases 
continuously without optimal control. In optimal 
control, the glucose concentration increases to around 
the 15‘ hour and then decreases continuously. Figure 


4b indicates insulin concentration before f{f) and 
after Ig{t) optimal control. The concentration before 


the optimal condition is continuous, while after the 
concentration increases, the concentration growth 
reduces from around the 15*" hour. The observation 
in Figure 4 shows that the optimal control, increasing 
the constant rate of insulin-independent glucose 
disappearance, effectively reduces the glucose 
concentration from the 15** hour. 


7. Conclusion 

The study presented a mathematical model 
formulating the concentration of glucose and insulin 
where extra glucose intake leading to extra 
stimulation of insulin production is considered. The 
study was motivated by the daily need to regulate 
glucose concentration in the body to manage Type-I 
diabetes. The well-posedness of the system indicated 
that the proposed model is applicable in real-life 
applications. Numerical simulation indicates that both 
glucose and insulin concentration continuously 
increases. This indicates a lack of optimality, a case 
that could persistently be type-I diabetes. A 
sensitivity analysis showed a constant rate of insulin- 
independent glucose disappearance is the most 
sensitive in the model. The study results observed that 
the extra glucose intake stimulates extra insulin 
production, causing the glucose concentration in the 
body to increase. This observation pointed to a need 


for optimal control based on the most sensitive 
parameter. We established optimal control based on a 
constant rate of insulin-independent glucose 
disappearance. The optimal control proved effective 
since the glucose curve decreases after 15** hour. 
Future studies need to consider including other 
parameters, such as glucose accelerated loss due to 
exercise in the system. 
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